home *** CD-ROM | disk | FTP | other *** search
/ Nebula 2 / Nebula Two.iso / Apps / Astro / ephem / Source / main.c < prev    next >
Encoding:
C/C++ Source or Header  |  1995-06-12  |  25.5 KB  |  1,040 lines
  1. /* main "ephem" program. 
  2.  * -------------------------------------------------------------------
  3.  * Copyright (c) 1990,1991,1992 by Elwood Charles Downey
  4.  * 
  5.  * Permission is granted to make and distribute copies of this program
  6.  * free of charge, provided the copyright notice and this permission
  7.  * notice are preserved on all copies.  All other rights reserved.
  8.  * -------------------------------------------------------------------
  9.  * set options.
  10.  * init screen and circumstances.
  11.  * enter infinite loop updating screen and allowing operator input.
  12.  */
  13.  
  14. #include <stdio.h>
  15. #include <ctype.h>
  16. #include <signal.h>
  17. #include <setjmp.h>
  18. #include <math.h>
  19. #ifdef VMS
  20. #include <stdlib.h>
  21. #endif
  22. #ifdef NeXT
  23. # include <stdlib.h>
  24. #endif
  25. #include "astro.h"
  26. #include "circum.h"
  27. #include "screen.h"
  28.  
  29. extern char *strncpy();
  30. extern char *getenv();
  31.  
  32. /* shorthands for fields of a Now structure, now.
  33.  * first undo the ones for a Now pointer from circum.h.
  34.  */
  35. #undef mjd
  36. #undef lat
  37. #undef lng
  38. #undef tz
  39. #undef temp
  40. #undef pressure
  41. #undef height
  42. #undef epoch
  43. #undef tznm
  44.  
  45. #define mjd    now.n_mjd
  46. #define lat    now.n_lat
  47. #define lng    now.n_lng
  48. #define tz    now.n_tz
  49. #define temp    now.n_temp
  50. #define pressure now.n_pressure
  51. #define height    now.n_height
  52. #define epoch    now.n_epoch
  53. #define tznm    now.n_tznm
  54.  
  55. static jmp_buf fpe_err_jmp;    /* used to recover from SIGFPE */
  56. static char *cfgfile;        /* !0 if -c used */
  57. static char cfgdef[] = "ephem.cfg"; /* default configuration file name */
  58. static Now now;        /* where when and how, right now */
  59. static double tminc;    /* hrs to inc time by each loop; RTC means use clock */
  60. static int nstep;    /* steps to go before stopping */
  61. static int spause;    /* secs to pause between steps */
  62. static int optwi;    /* set when want to display dawn/dusk/len-of-night */
  63. static int oppl;    /* mask of (1<<planet) bits; set when want to show it */
  64.  
  65. main (ac, av)
  66. int ac;
  67. char *av[];
  68. {
  69.     void bye();
  70.     void on_fpe();
  71.     static char freerun[] =
  72.         "Running... press any key to stop to make changes.";
  73.     static char prmpt[] =
  74. "Move to another field, RETURN to change this field, ? for help, or q to run";
  75.     static char hlp[] =
  76.     "arrow keys move to field; any key stops running; ^d exits; ^l redraws";
  77.     int fld = rcfpack(R_NSTEP, C_NSTEPV, 0); /* initial cursor loc */
  78.     int sflag = 0;    /* not silent, by default */
  79.     int one = 1;    /* use a variable so optimizer doesn't get disabled */
  80.     int srchdone = 0; /* true when search funcs say so */
  81.     int newcir = 2;    /* set when circumstances change - means don't tminc */
  82.  
  83.     while ((--ac > 0) && (**++av == '-')) {
  84.         char *s;
  85.         for (s = *av+1; *s != '\0'; s++)
  86.         switch (*s) {
  87.         case 's': /* no credits "silent" (don't publish this) */
  88.             sflag++;
  89.             break;
  90.         case 'c': /* set name of config file to use */
  91.             if (--ac <= 0) usage("-c but no config file");
  92.             cfgfile = *++av;
  93.             break;
  94.         case 'd': /* set alternate database file name */
  95.             if (--ac <= 0) usage("-d but no database file");
  96.             obj_setdbfilename (*++av);
  97.             break;
  98.         default:
  99.             usage("Bad - option");
  100.         }
  101.     }
  102.  
  103.     if (!sflag)
  104.         credits();
  105.  
  106.     /* fresh screen.
  107.      * crack config file, THEN args so args may override.
  108.      */
  109.     c_erase();
  110.     read_cfgfile ();
  111.     read_fieldargs (ac, av);
  112.  
  113.     /* set up to clean up screen and tty if interrupted.
  114.      * also set up to stop if get floating error.
  115.      */
  116.     (void) signal (SIGINT, bye);
  117.     (void) signal (SIGFPE, on_fpe);
  118.  
  119.     /* update screen forever (until QUIT) */
  120.     while (one) {
  121.  
  122.         /* if get a floating error, longjmp() here and stop looping */
  123.         if (setjmp (fpe_err_jmp))
  124.         nstep = 0;
  125.         else {
  126.         nstep -= 1;
  127.  
  128.         /* recalculate everything and update all the fields */
  129.         redraw_screen (newcir);
  130.         mm_newcir (0);
  131.  
  132.         /* let searching functions change tminc and check for done */
  133.         srchdone = srch_eval (mjd, &tminc) < 0;
  134.         print_tminc(0);    /* to show possibly new search increment */
  135.  
  136.         /* update plot and listing files, now that all fields are up
  137.          * to date and search function has been evaluated.
  138.          */
  139.         plot();
  140.         listing();
  141.  
  142.         /* handle spause if we are really looping */
  143.         if (nstep > 0)
  144.             slp_sync();
  145.         }
  146.  
  147.         /* stop loop to allow op to change parameters:
  148.          * if a search evaluation converges (or errors out),
  149.          * or if steps are done,
  150.          * or if op hits any key.
  151.          */
  152.         newcir = 0;
  153.         if (srchdone || nstep <= 0 || (chk_char()==0 && read_char()!=0)) {
  154.         int nfld;
  155.  
  156.         /* update screen with the current stuff if stopped during
  157.          * unattended plotting or listing since last redraw_screen()
  158.          * didn't.
  159.          */
  160.         if ((plot_ison() || listing_ison()) && nstep > 0)
  161.             redraw_screen (1);
  162.  
  163.         /* return nstep to default of 1 */
  164.         if (nstep <= 0) {
  165.             nstep = 1;
  166.             print_nstep (0);
  167.         }
  168.  
  169.         /* change fields until END.
  170.          * update all time fields if any are changed
  171.          * and print NEW CIRCUMSTANCES if any have changed.
  172.          * QUIT causes bye() to be called and we never return.
  173.          */
  174.         while(nfld = sel_fld(fld,alt_menumask()|F_CHG,prmpt,hlp)) {
  175.             if (chg_fld ((char *)0, &nfld)) {
  176.             mm_now (&now, 1);
  177.             mm_newcir(1);
  178.             newcir = 1;
  179.             }
  180.             fld = nfld;
  181.         }
  182.         if (nstep > 1)
  183.             f_prompt (freerun);
  184.         }
  185.  
  186.         /* increment time only if op didn't change cirumstances */
  187.         if (!newcir)
  188.         inc_mjd (&now, tminc);
  189.     }
  190.  
  191.     return (0);
  192. }
  193.  
  194. /* read in ephem's configuration file, if any.
  195.  * if errors in file, call usage() (which exits).
  196.  * if use -d, require it; else try $EPHEMCFG and ephem.cfg but don't
  197.  *   complain if can't find these since, after all, one is not required.
  198.  * skip all lines that doesn't begin with an alpha char.
  199.  */
  200. static
  201. read_cfgfile()
  202. {
  203.     char buf[128];
  204.     FILE *fp;
  205.     char *fn;
  206.  
  207.     /* open the config file. 
  208.      * only REQUIRED if used -d option.
  209.      * if succcessful, fn points to file name.
  210.      */
  211.     if (cfgfile) {
  212.         fn = cfgfile;
  213.         fp = fopen (fn, "r");
  214.         if (!fp) {
  215.         (void) sprintf (buf, "Can not open %s", fn);
  216.         usage (buf);    /* does not return */
  217.         }
  218.     } else {
  219.         fn = getenv ("EPHEMCFG");
  220.         if (!fn)
  221.         fn = cfgdef;
  222.     }
  223.     fp = fopen (fn, "r");
  224.     if (!fp)
  225.         return;    /* oh well; after all, it's not required */
  226.  
  227.     while (fgets (buf, sizeof(buf), fp)) {
  228.         if (!isalpha(buf[0]))
  229.         continue;
  230.         buf[strlen(buf)-1] = '\0';        /* discard trailing \n */
  231.         if (crack_fieldset (buf) < 0) {
  232.         char why[NC];
  233.         (void) sprintf (why, "Bad field spec in %s: %s\n", fn, buf);
  234.         usage (why);
  235.         }
  236.     }
  237.     (void) fclose (fp);
  238. }
  239.  
  240.  
  241. /* draw all the stuff on the screen, using the current menu.
  242.  * if how_much == 0 then just update fields that need it;
  243.  * if how_much == 1 then redraw all fields;
  244.  * if how_much == 2 then erase the screen and redraw EVERYTHING.
  245.  */
  246. redraw_screen (how_much)
  247. int how_much;
  248. {
  249.     if (how_much == 2)
  250.         c_erase();
  251.  
  252.     /* print the single-step message if this is the last loop */
  253.     if (nstep < 1)
  254.         print_updating();
  255.  
  256.     if (how_much == 2) {
  257.         mm_borders();
  258.         mm_labels();
  259.         srch_prstate(1);
  260.         plot_prstate(1);
  261.         listing_prstate(1);
  262.     }
  263.  
  264.     /* if just updating changed fields while plotting or listing
  265.      * unattended then suppress most screen updates except
  266.      * always show nstep to show plot loops to go and
  267.      * always show tminc to show search convergence progress.
  268.      */
  269.     print_nstep(how_much);
  270.     print_tminc(how_much);
  271.     print_spause(how_much);
  272.     if (how_much == 0 && (plot_ison() || listing_ison()) && nstep > 0)
  273.         f_off();
  274.  
  275.     /* print all the time-related fields */
  276.     mm_now (&now, how_much);
  277.  
  278.     if (optwi)
  279.         mm_twilight (&now, how_much);
  280.  
  281.     /* print stuff on bottom menu */
  282.     print_alt (how_much);
  283.  
  284.     f_on();
  285. }
  286.  
  287. /* clean up and exit.
  288.  */
  289. void
  290. bye()
  291. {
  292.     c_erase();
  293.     byetty();
  294.     exit (0);
  295. }
  296.  
  297. /* this gets called when a floating point error occurs.
  298.  * we force a jump back into main() with looping terminated.
  299.  */
  300. static
  301. void
  302. on_fpe()
  303. {
  304.     extern void longjmp();
  305.  
  306.     (void) signal (SIGFPE, on_fpe);
  307.     f_msg ("Floating point error has occurred - computations aborted.");
  308.     longjmp (fpe_err_jmp, 1);
  309. }
  310.  
  311. usage(why)
  312. char *why;
  313. {
  314.     /* don't advertise -s (silent) option */
  315.     c_erase();
  316.     f_string (1, 1, why);
  317.     f_string (2, 1,
  318.         "usage: [-c <configfile>] [-d <database>] [field=value ...]\r\n");
  319.     byetty();
  320.     exit (1);
  321. } 
  322.  
  323. /* Process the field specs from the command line.
  324.  * if trouble call usage() (which exits).
  325.  */
  326. static
  327. read_fieldargs (ac, av)
  328. int ac;        /* number of such specs */
  329. char *av[];    /* array of strings in form <field_name value> */
  330. {
  331.     while (--ac >= 0) {
  332.         char *fs = *av++;
  333.         if (crack_fieldset (fs) < 0) {
  334.         char why[NC];
  335.         (void) sprintf (why, "Bad command line spec: %.*s",
  336.                             sizeof(why)-26, fs);
  337.         usage (why);
  338.         }
  339.     }
  340. }
  341.  
  342. /* process a field spec in buf, either from config file or argv.
  343.  * return 0 if recognized ok, else -1.
  344.  */
  345. static
  346. crack_fieldset (buf)
  347. char *buf;
  348. {
  349. #define    ARRAY_SIZ(a)    (sizeof(a)/sizeof((a)[0]))
  350. #define    MAXKW        6    /* longest keyword, not counting trailing 0 */
  351.     /* N.B. index of item is its case value, below.
  352.      * N.B. if add an item, keep it no longer than MAXKW chars.
  353.      */
  354.     static char keywords[][MAXKW+1] = {
  355.         /*  0 */    "LAT",
  356.         /*  1 */    "LONG",
  357.         /*  2 */    "UT",
  358.         /*  3 */    "UD",
  359.         /*  4 */    "TZONE",
  360.         /*  5 */    "TZNAME",
  361.         /*  6 */    "HEIGHT",
  362.         /*  7 */    "NSTEP",
  363.         /*  8 */    "PAUSE",
  364.         /*  9 */    "STPSZ",
  365.         /* 10 */    "TEMP",
  366.         /* 11 */    "PRES",
  367.         /* 12 */    "EPOCH",
  368.         /* 13 */    "JD",
  369.         /* 14 */    "OBJX",
  370.         /* 15 */    "OBJY",
  371.         /* 16 */    "PROPTS",
  372.         /* 17 */    "MENU"
  373.     };
  374.     int i;
  375.     int l;
  376.     int f;
  377.  
  378.     for (i = 0; i < ARRAY_SIZ(keywords); i++)
  379.         if (strncmp (keywords[i], buf, l = strlen(keywords[i])) == 0) {
  380.         buf += l+1;    /* skip keyword and its subsequent delimiter */
  381.         break;
  382.         }
  383.  
  384.     switch (i) {
  385.     case 0: f = rcfpack (R_LAT,C_LATV,0); (void) chg_fld (buf, &f);
  386.         break;
  387.     case 1: f = rcfpack (R_LONG,C_LONGV,0), (void) chg_fld (buf, &f);
  388.         break;
  389.     case 2: f = rcfpack (R_UT,C_UTV,0), (void) chg_fld (buf, &f);
  390.         break;
  391.     case 3: f = rcfpack (R_UD,C_UD,0), (void) chg_fld (buf, &f);
  392.         break;
  393.     case 4: f = rcfpack (R_TZONE,C_TZONEV,0), (void) chg_fld (buf, &f);
  394.         break;
  395.     case 5: f = rcfpack (R_TZN,C_TZN,0), (void) chg_fld (buf, &f);
  396.         break;
  397.     case 6: f = rcfpack (R_HEIGHT,C_HEIGHTV,0), (void) chg_fld (buf, &f);
  398.         break;
  399.     case 7: f = rcfpack (R_NSTEP,C_NSTEPV,0), (void) chg_fld (buf, &f);
  400.         break;
  401.     case 8: f = rcfpack (R_PAUSE,C_PAUSEV,0), (void) chg_fld (buf, &f);
  402.         break;
  403.     case 9: f = rcfpack (R_STPSZ,C_STPSZV,0), (void) chg_fld (buf, &f);
  404.         break;
  405.     case 10: f = rcfpack (R_TEMP,C_TEMPV,0), (void) chg_fld (buf, &f);
  406.         break;
  407.     case 11: f = rcfpack (R_PRES,C_PRESV,0), (void) chg_fld (buf, &f);
  408.         break;
  409.     case 12: f = rcfpack (R_EPOCH,C_EPOCHV,0), (void) chg_fld (buf, &f);
  410.         break;
  411.     case 13: f = rcfpack (R_JD,C_JDV,0), (void) chg_fld (buf, &f);
  412.         break;
  413.     case 14: (void) obj_filelookup (OBJX, buf);
  414.         break;
  415.     case 15: (void) obj_filelookup (OBJY, buf);
  416.         break;
  417.     case 16:
  418.         if (buf[-1] != '+')
  419.         optwi = oppl = 0;
  420.         while (*buf)
  421.         switch (*buf++) {
  422.         case 'T': optwi = 1; break;
  423.         case 'S': oppl |= (1<<SUN); break;
  424.         case 'M': oppl |= (1<<MOON); break;
  425.         case 'e': oppl |= (1<<MERCURY); break;
  426.         case 'v': oppl |= (1<<VENUS); break;
  427.         case 'm': oppl |= (1<<MARS); break;
  428.         case 'j': case 'J': oppl |= (1<<JUPITER); break;
  429.         case 's': oppl |= (1<<SATURN); break;
  430.         case 'u': oppl |= (1<<URANUS); break;
  431.         case 'n': oppl |= (1<<NEPTUNE); break;
  432.         case 'p': oppl |= (1<<PLUTO); break;
  433.         case 'x': oppl |= (1<<OBJX); obj_on(OBJX); break;
  434.         case 'y': oppl |= (1<<OBJY); obj_on(OBJY); break;
  435.         }
  436.         break;
  437.     case 17:
  438.         if (strncmp (buf, "DATA", 4) == 0)
  439.         altmenu_init (F_MNU1);
  440.         else if (strncmp (buf, "RISET", 5) == 0)
  441.         altmenu_init (F_MNU2);
  442.         else if (strncmp (buf, "SEP", 3) == 0)
  443.         altmenu_init (F_MNU3);
  444.         else if (strncmp (buf, "JUP", 3) == 0)
  445.         altmenu_init (F_MNUJ);
  446.         break;
  447.     default:
  448.         return (-1);
  449.     }
  450.     return (0);
  451. }
  452.  
  453. /* react to the field at *fld according to the optional string input at bp.
  454.  * if bp is != 0 use it, else issue read_line() and use buffer.
  455.  * then sscanf the buffer and update the corresponding (global) variable(s)
  456.  * or do whatever a pick at that field should do.
  457.  * we might also change *fld if we want to change the current cursor location.
  458.  * return 1 if we change a field that invalidates any of the times or
  459.  * to update all related fields.
  460.  */
  461. static
  462. chg_fld (bp, fld)
  463. char *bp;
  464. int *fld;
  465. {
  466.     char buf[NC];
  467.     int deghrs = 0, mins = 0, secs = 0;
  468.     int new = 0;
  469.  
  470.     /* switch on just the row/col portion */
  471.     switch (unpackrc(*fld)) {
  472.     case rcfpack (R_ALTM, C_ALTM, 0):
  473.         if (altmenu_setup() == 0) {
  474.         print_updating();
  475.         alt_erase();
  476.         print_alt(2);
  477.         }
  478.         break;
  479.     case rcfpack (R_JD, C_JDV, 0):
  480.         if (!bp) {
  481.         static char p[] = "Julian Date (or n for Now): ";
  482.         f_prompt (p);
  483.         if (read_line (buf, PW-sizeof(p)) <= 0)
  484.             break;
  485.         bp = buf;
  486.         }
  487.         if (bp[0] == 'n' || bp[0] == 'N')
  488.         time_fromsys (&now);
  489.         else
  490.         mjd = atof(bp) - 2415020L;
  491.         set_t0 (&now);
  492.         new = 1;
  493.         break;
  494.     case rcfpack (R_UD, C_UD, 0):
  495.         if (!bp) {
  496.         static char p[] = "utc date (m/d/y, or year.d, or n for Now): ";
  497.         f_prompt (p);
  498.         if (read_line (buf, PW-sizeof(p)) <= 0)
  499.             break;
  500.         bp = buf;
  501.         }
  502.         if (bp[0] == 'n' || bp[0] == 'N')
  503.         time_fromsys (&now);
  504.         else {
  505.         if (decimal_year(bp)) {
  506.             double y = atof (bp);
  507.             year_mjd (y, &mjd);
  508.         } else {
  509.             double day, newmjd0;
  510.             int month, year;
  511.             mjd_cal (mjd, &month, &day, &year); /* init with now */
  512.             f_sscandate (bp, &month, &day, &year);
  513.             cal_mjd (month, day, year, &newmjd0);
  514.             /* if don't give a fractional part to days
  515.              * then retain current hours.
  516.              */
  517.             if ((long)day == day)
  518.             mjd = newmjd0 + mjd_hr(mjd)/24.0;
  519.             else
  520.             mjd = newmjd0;
  521.         }
  522.         }
  523.         set_t0 (&now);
  524.         new = 1;
  525.         break;
  526.     case rcfpack (R_UT, C_UTV, 0):
  527.         if (!bp) {
  528.         static char p[] = "utc time (h:m:s, or n for Now): ";
  529.         f_prompt (p);
  530.         if (read_line (buf, PW-sizeof(p)) <= 0)
  531.             break;
  532.         bp = buf;
  533.         }
  534.         if (bp[0] == 'n' || bp[0] == 'N')
  535.         time_fromsys (&now);
  536.         else {
  537.         double newutc = (mjd-mjd_day(mjd)) * 24.0;
  538.         f_dec_sexsign (newutc, °hrs, &mins, &secs);
  539.         f_sscansex (bp, °hrs, &mins, &secs);
  540.         sex_dec (deghrs, mins, secs, &newutc);
  541.         mjd = mjd_day(mjd) + newutc/24.0;
  542.         }
  543.         set_t0 (&now);
  544.         new = 1;
  545.         break;
  546.     case rcfpack (R_LD, C_LD, 0):
  547.         if (!bp) {
  548.         static char p[] = "local date (m/d/y, or year.d, n for Now): ";
  549.         f_prompt (p);
  550.         if (read_line (buf, PW-sizeof(p)) <= 0)
  551.             break;
  552.         bp = buf;
  553.         }
  554.         if (bp[0] == 'n' || bp[0] == 'N')
  555.         time_fromsys (&now);
  556.         else {
  557.         if (decimal_year(bp)) {
  558.             double y = atof (bp);
  559.             year_mjd (y, &mjd);
  560.             mjd += tz/24.0;
  561.         } else {
  562.             double day, newlmjd0;
  563.             int month, year;
  564.             mjd_cal (mjd-tz/24.0, &month, &day, &year); /* now */
  565.             f_sscandate (bp, &month, &day, &year);
  566.             cal_mjd (month, day, year, &newlmjd0);
  567.             /* if don't give a fractional part to days
  568.              * then retain current hours.
  569.              */
  570.             if ((long)day == day)
  571.             mjd = newlmjd0 + mjd_hr(mjd-tz/24.0)/24.0;
  572.             else
  573.             mjd = newlmjd0;
  574.             mjd += tz/24.0;
  575.         }
  576.         }
  577.         set_t0 (&now);
  578.         new = 1;
  579.         break;
  580.     case rcfpack (R_LT, C_LT, 0):
  581.         if (!bp) {
  582.         static char p[] = "local time (h:m:s, or n for Now): ";
  583.         f_prompt (p);
  584.         if (read_line (buf, PW-sizeof(p)) <= 0)
  585.             break;
  586.         bp = buf;
  587.         }
  588.         if (bp[0] == 'n' || bp[0] == 'N')
  589.         time_fromsys (&now);
  590.         else {
  591.         double newlt = (mjd-mjd_day(mjd)) * 24.0 - tz;
  592.         range (&newlt, 24.0);
  593.         f_dec_sexsign (newlt, °hrs, &mins, &secs);
  594.         f_sscansex (bp, °hrs, &mins, &secs);
  595.         sex_dec (deghrs, mins, secs, &newlt);
  596.         mjd = mjd_day(mjd-tz/24.0) + (newlt + tz)/24.0;
  597.         }
  598.         set_t0 (&now);
  599.         new = 1;
  600.         break;
  601.     case rcfpack (R_LST, C_LSTV, 0):
  602.         if (!bp) {
  603.         static char p[] = "local sidereal time (h:m:s, or n for Now): ";
  604.         f_prompt (p);
  605.         if (read_line (buf, PW-sizeof(p)) <= 0)
  606.             break;
  607.         bp = buf;
  608.         }
  609.         if (bp[0] == 'n' || bp[0] == 'N')
  610.         time_fromsys (&now);
  611.         else {
  612.         double lst, utc;
  613.         now_lst (&now, &lst);
  614.         f_dec_sexsign (lst, °hrs, &mins, &secs);
  615.         f_sscansex (bp, °hrs, &mins, &secs);
  616.         sex_dec (deghrs, mins, secs, &lst);
  617.         lst -= radhr(lng); /* convert to gst */
  618.         range (&lst, 24.0);
  619.         gst_utc (mjd_day(mjd), lst, &utc);
  620.         mjd = mjd_day(mjd) + utc/24.0;
  621.         }
  622.         set_t0 (&now);
  623.         new = 1;
  624.         break;
  625.     case rcfpack (R_TZN, C_TZN, 0):
  626.         if (!bp) {
  627.         static char p[] = "timezone abbreviation (3 char max): ";
  628.         f_prompt (p);
  629.         if (read_line (buf, 3) <= 0)
  630.             break;
  631.         bp = buf;
  632.         }
  633.         (void) strncpy (tznm, bp, sizeof(tznm)-1);
  634.         new = 1;
  635.         break;
  636.     case rcfpack (R_TZONE, C_TZONEV, 0):
  637.         if (!bp) {
  638.         static char p[] = "hours behind utc: ";
  639.         f_prompt (p);
  640.         if (read_line (buf, PW-sizeof(p)) <= 0)
  641.             break;
  642.         bp = buf;
  643.         }
  644.         f_dec_sexsign (tz, °hrs, &mins, &secs);
  645.         f_sscansex (bp, °hrs, &mins, &secs);
  646.         sex_dec (deghrs, mins, secs, &tz);
  647.         new = 1;
  648.         break;
  649.     case rcfpack (R_LONG, C_LONGV, 0):
  650.         if (!bp) {
  651.         static char p[] = "longitude (+ west) (d:m:s): ";
  652.         f_prompt (p);
  653.         if (read_line (buf, PW-sizeof(p)) <= 0)
  654.             break;
  655.         bp = buf;
  656.         }
  657.         f_dec_sexsign (-raddeg(lng), °hrs, &mins, &secs);
  658.         f_sscansex (bp, °hrs, &mins, &secs);
  659.         sex_dec (deghrs, mins, secs, &lng);
  660.         lng = degrad (-lng);         /* want - radians west */
  661.         new = 1;
  662.         break;
  663.     case rcfpack (R_LAT, C_LATV, 0):
  664.         if (!bp) {
  665.         static char p[] = "latitude (+ north) (d:m:s): ";
  666.         f_prompt (p);
  667.         if (read_line (buf, PW-sizeof(p)) <= 0)
  668.             break;
  669.         bp = buf;
  670.         }
  671.         f_dec_sexsign (raddeg(lat), °hrs, &mins, &secs);
  672.         f_sscansex (bp, °hrs, &mins, &secs);
  673.         sex_dec (deghrs, mins, secs, &lat);
  674.         lat = degrad (lat);
  675.         new = 1;
  676.         break;
  677.     case rcfpack (R_HEIGHT, C_HEIGHTV, 0):
  678.         if (!bp) {
  679.         static char p[] = "height above sea level (ft): ";
  680.         f_prompt (p);
  681.         if (read_line (buf, PW-sizeof(p)) <= 0)
  682.             break;
  683.         bp = buf;
  684.         }
  685.         if (sscanf (bp, "%lf", &height) == 1) {
  686.         height /= 2.093e7; /*convert ft to earth radii above sea level*/
  687.         new = 1;
  688.         }
  689.         break;
  690.     case rcfpack (R_NSTEP, C_NSTEPV, 0):
  691.         if (!bp) {
  692.         static char p[] = "number of steps to run: ";
  693.         f_prompt (p);
  694.         if (read_line (buf, 8) <= 0)
  695.             break;
  696.         bp = buf;
  697.         }
  698.         (void) sscanf (bp, "%d", &nstep);
  699.         print_nstep (0);
  700.         break;
  701.     case rcfpack (R_PAUSE, C_PAUSEV, 0):
  702.         if (!bp) {
  703.         static char p[] = "seconds to pause between steps: ";
  704.         f_prompt (p);
  705.         if (read_line (buf, 8) <= 0)
  706.             break;
  707.         bp = buf;
  708.         }
  709.         (void) sscanf (bp, "%d", &spause);
  710.         print_spause (0);
  711.         break;
  712.     case rcfpack (R_TEMP, C_TEMPV, 0):
  713.         if (!bp) {
  714.         static char p[] = "temperature (deg.F): ";
  715.         f_prompt (p);
  716.         if (read_line (buf, PW-sizeof(p)) <= 0)
  717.             break;
  718.         bp = buf;
  719.         }
  720.         if (sscanf (bp, "%lf", &temp) == 1) {
  721.         temp = 5./9.*(temp - 32.0);    /* want degs C */
  722.         new = 1;
  723.         }
  724.         break;
  725.     case rcfpack (R_PRES, C_PRESV, 0):
  726.         if (!bp) {
  727.         static char p[] =
  728.             "atmos pressure (in. Hg; 0 for no refraction correction): ";
  729.         f_prompt (p);
  730.         if (read_line (buf, PW-sizeof(p)) <= 0)
  731.             break;
  732.         bp = buf;
  733.         }
  734.         if (sscanf (bp, "%lf", &pressure) == 1) {
  735.         pressure *= 33.86;        /* want mBar */
  736.         new = 1;
  737.         }
  738.         break;
  739.     case rcfpack (R_EPOCH, C_EPOCHV, 0):
  740.         if (!bp) {
  741.         static char p[] = "epoch (year, or e for Equinox of Date): ";
  742.         f_prompt (p);
  743.         if (read_line (buf, PW-strlen(p)) <= 0)
  744.             break;
  745.         bp = buf;
  746.         }
  747.         if (bp[0] == 'e' || bp[0] == 'E')
  748.         epoch = EOD;
  749.         else {
  750.         double e;
  751.         e = atof(bp);
  752.         year_mjd (e, &epoch);
  753.         }
  754.         new = 1;
  755.         break;
  756.     case rcfpack (R_STPSZ, C_STPSZV, 0):
  757.         if (!bp) {
  758.         static char p[] =
  759.             "step size increment (h:m:s, or <x>d for x days, or r for RTC): ";
  760.         f_prompt (p);
  761.         if (read_line (buf, PW-sizeof(p)) <= 0)
  762.             break;
  763.         bp = buf;
  764.         }
  765.         if (bp[0] == 'r' || bp[0] == 'R')
  766.         tminc = RTC;
  767.         else {
  768.         int last = strlen (bp) - 1;
  769.         if (bp[last] == 'd') {
  770.             /* ends in d so treat as a number of days */
  771.             double x;
  772.             if (sscanf (bp, "%lf", &x) == 1)
  773.             tminc = x * 24.0;
  774.         } else {
  775.             if (tminc == RTC)
  776.             deghrs = mins = secs = 0;
  777.             else
  778.             f_dec_sexsign (tminc, °hrs, &mins, &secs);
  779.             f_sscansex (bp, °hrs, &mins, &secs);
  780.             sex_dec (deghrs, mins, secs, &tminc);
  781.         }
  782.         }
  783.         print_tminc(0);
  784.         set_t0 (&now);
  785.         break;
  786.     case rcfpack (R_PLOT, C_PLOT, 0):
  787.         plot_setup();
  788.         if (plot_ison())
  789.         new = 1;
  790.         break;
  791.     case rcfpack (R_LISTING, C_LISTING, 0):
  792.         listing_setup();
  793.         if (listing_ison())
  794.         new = 1;
  795.         break;
  796.     case rcfpack (R_WATCH, C_WATCH, 0):
  797.         watch (&now, tminc, oppl);
  798.         /* set new reference time to what watch left it.
  799.          * no need to set new since watch just did a redraw.
  800.          */
  801.         set_t0 (&now);
  802.         break;
  803.     case rcfpack (R_DAWN, C_DAWN, 0):
  804.     case rcfpack (R_DUSK, C_DUSK, 0):
  805.     case rcfpack (R_LON, C_LON, 0):
  806.         if (optwi ^= 1) {
  807.         print_updating();
  808.         mm_twilight (&now, 1);
  809.         } else {
  810.         f_blanks (R_DAWN, C_DAWNV, 5);
  811.         f_blanks (R_DUSK, C_DUSKV, 5);
  812.         f_blanks (R_LON, C_LONV, 5);
  813.         }
  814.         break;
  815.     case rcfpack (R_SRCH, C_SRCH, 0):
  816.         srch_setup();
  817.         if (srch_ison())
  818.         new = 1;
  819.         break;
  820.     case rcfpack (R_SUN, C_OBJ, 0):
  821.         toggle_body (SUN);
  822.         break;
  823.     case rcfpack (R_SUN, C_CONSTEL, 0):
  824.         if (oppl & (1<<SUN))
  825.         constellation_msg (SUN, &now);
  826.         break;
  827.     case rcfpack (R_MOON, C_OBJ, 0):
  828.         toggle_body (MOON);
  829.         break;
  830.     case rcfpack (R_MOON, C_CONSTEL, 0):
  831.         if (oppl & (1<<MOON))
  832.         constellation_msg (MOON, &now);
  833.         break;
  834.     case rcfpack (R_MERCURY, C_OBJ, 0):
  835.         toggle_body (MERCURY);
  836.         break;
  837.     case rcfpack (R_MERCURY, C_CONSTEL, 0):
  838.         if (oppl & (1<<MERCURY))
  839.         constellation_msg (MERCURY, &now);
  840.         break;
  841.     case rcfpack (R_VENUS, C_OBJ, 0):
  842.         toggle_body (VENUS);
  843.         break;
  844.     case rcfpack (R_VENUS, C_CONSTEL, 0):
  845.         if (oppl & (1<<VENUS))
  846.         constellation_msg (VENUS, &now);
  847.         break;
  848.     case rcfpack (R_MARS, C_OBJ, 0):
  849.         toggle_body (MARS);
  850.         break;
  851.     case rcfpack (R_MARS, C_CONSTEL, 0):
  852.         if (oppl & (1<<MARS))
  853.         constellation_msg (MARS, &now);
  854.         break;
  855.     case rcfpack (R_JUPITER, C_OBJ, 0):
  856.         toggle_body (JUPITER);
  857.         break;
  858.     case rcfpack (R_JUPITER, C_CONSTEL, 0):
  859.         if (oppl & (1<<JUPITER))
  860.         constellation_msg (JUPITER, &now);
  861.         break;
  862.     case rcfpack (R_JUPITER, C_XTRA, 0):
  863.         if (oppl & (1<<JUPITER)) {
  864.         print_updating();
  865.         alt_erase();
  866.         altmenu_init (F_MNUJ);
  867.         print_alt (2);
  868.         *fld = rcfpack(R_NSTEP, C_NSTEPV, 0);
  869.         }
  870.         break;
  871.     case rcfpack (R_SATURN, C_OBJ, 0):
  872.         toggle_body (SATURN);
  873.         break;
  874.     case rcfpack (R_SATURN, C_CONSTEL, 0):
  875.         if (oppl & (1<<SATURN))
  876.         constellation_msg (SATURN, &now);
  877.         break;
  878.     case rcfpack (R_URANUS, C_OBJ, 0):
  879.         toggle_body (URANUS);
  880.         break;
  881.     case rcfpack (R_URANUS, C_CONSTEL, 0):
  882.         if (oppl & (1<<URANUS))
  883.         constellation_msg (URANUS, &now);
  884.         break;
  885.     case rcfpack (R_NEPTUNE, C_OBJ, 0):
  886.         toggle_body (NEPTUNE);
  887.         break;
  888.     case rcfpack (R_NEPTUNE, C_CONSTEL, 0):
  889.         if (oppl & (1<<NEPTUNE))
  890.         constellation_msg (NEPTUNE, &now);
  891.         break;
  892.     case rcfpack (R_PLUTO, C_OBJ, 0):
  893.         toggle_body (PLUTO);
  894.         break;
  895.     case rcfpack (R_PLUTO, C_CONSTEL, 0):
  896.         if (oppl & (1<<PLUTO))
  897.         constellation_msg (PLUTO, &now);
  898.         break;
  899.     case rcfpack (R_OBJX, C_OBJ, 0):
  900.         /* this might change which columns are used so erase all when
  901.          * returns and redraw if still on.
  902.          */
  903.         obj_setup (OBJX);
  904.         alt_nobody (OBJX);
  905.         if (obj_ison (OBJX)) {
  906.         oppl |= 1 << OBJX;
  907.         print_updating();
  908.         alt_body (OBJX, 1, &now);
  909.         } else
  910.         oppl &= ~(1 << OBJX);    /* already erased; just clear flag */
  911.         break;
  912.     case rcfpack (R_OBJX, C_CONSTEL, 0):
  913.         if (oppl & (1<<OBJX))
  914.         constellation_msg (OBJX, &now);
  915.         break;
  916.     case rcfpack (R_OBJY, C_OBJ, 0):
  917.         /* this might change which columns are used so erase all when
  918.          * returns and redraw if still on.
  919.          */
  920.         obj_setup (OBJY);
  921.         alt_nobody (OBJY);
  922.         if (obj_ison (OBJY)) {
  923.         oppl |= 1 << OBJY;
  924.         print_updating();
  925.         alt_body (OBJY, 1, &now);
  926.         } else
  927.         oppl &= ~(1 << OBJY);    /* already erased; just clear flag */
  928.         break;
  929.     case rcfpack (R_OBJY, C_CONSTEL, 0):
  930.         if (oppl & (1<<OBJY))
  931.         constellation_msg (OBJY, &now);
  932.         break;
  933.     }
  934.  
  935.     return (new);
  936. }
  937.  
  938. static
  939. print_tminc(force)
  940. int force;
  941. {
  942.     static double last = -123.456;    /* anything unlikely */
  943.  
  944.     if (force || tminc != last) {
  945.         if (tminc == RTC)
  946.         f_string (R_STPSZ, C_STPSZV, " RT CLOCK");
  947.         else if (fabs(tminc) >= 24.0)
  948.         f_double (R_STPSZ, C_STPSZV, "%6.4g dy", tminc/24.0);
  949.         else
  950.         f_signtime (R_STPSZ, C_STPSZV, tminc);
  951.         last = tminc;
  952.     }
  953. }
  954.  
  955. /* print stuff on bottom menu */
  956. static
  957. print_alt (howmuch)
  958. int howmuch;
  959. {
  960.     if (howmuch == 2)
  961.         alt_labels();
  962.     if (alt_menumask() == F_MNUJ)
  963.         altj_display (howmuch, &now);
  964.     else {
  965.         int p;
  966.         for (p = nxtbody(-1); p != -1; p = nxtbody(p))
  967.         if (oppl & (1<<p))
  968.             alt_body (p, howmuch, &now);
  969.     }
  970. }
  971.  
  972. print_updating()
  973. {
  974.     f_prompt ("Updating...");
  975. }
  976.  
  977. static
  978. print_nstep(force)
  979. int force;
  980. {
  981.     static int last;
  982.  
  983.     if (force || nstep != last) {
  984.         char buf[16];
  985.         (void) sprintf (buf, "%8d", nstep);
  986.         f_string (R_NSTEP, C_NSTEPV, buf);
  987.         last = nstep;
  988.     }
  989. }
  990.  
  991. static
  992. print_spause(force)
  993. int force;
  994. {
  995.     static int last;
  996.  
  997.     if (force || spause != last) {
  998.         char buf[16];
  999.         (void) sprintf (buf, "%8d", spause);
  1000.         f_string (R_PAUSE, C_PAUSEV, buf);
  1001.         last = spause;
  1002.     }
  1003. }
  1004.  
  1005. /* if not plotting/listing/searching then sleep spause seconds.
  1006.  * if time is being based on the real-time clock, sync on the next
  1007.  *   integral multiple of spause seconds after the minute.
  1008.  * check for keyboard action once each second to let it break out early.
  1009.  */
  1010. slp_sync()
  1011. {
  1012.     extern long time();
  1013.  
  1014.     if (spause > 0 && !plot_ison() && !srch_ison() && !listing_ison()) {
  1015.         int n;
  1016.         if (tminc == RTC) {
  1017.         long t;
  1018.         (void) time (&t);
  1019.         n = spause - (t % spause);
  1020.         } else 
  1021.         n = spause;
  1022.         while (--n >= 0)
  1023.         if (chk_char() == 0)
  1024.             break;
  1025.         else
  1026.             (void) sleep (1);
  1027.     }
  1028. }
  1029.  
  1030. static
  1031. toggle_body (p)
  1032. int p;
  1033. {
  1034.     if ((oppl ^= (1<<p)) & (1<<p)) {
  1035.         print_updating();
  1036.         alt_body (p, 1, &now);
  1037.     } else
  1038.         alt_nobody (p);
  1039. }
  1040.